Accepted for publication in Math. Proc. Camb. Phil. Soc. 



1 



Complex asymptotics of Poincare functions and properties of 

Julia sets 

By GREGORY DERFEL 
Department of Mathematics and Computer Science, 
Ben Gurion University of the Negev, Beer Sheva 84-105, Israel 
e-mail: derfel@math.bgu.ac.il 

PETER J. GRABNERf 
Institut fur Analysis und Computational Number Theory (Math A), 
Technische Universitat Graz, Steyrergasse 30, 8010 Graz, Austria 
e-mail: peter .grabnerOtugraz. at 

and FRITZ VOGL 

Institut fur Analysis und Scientific Computing, Technische Universitat Wien, 
Wiedner Hauptstrafle 8-10, 1040 Wien, Austria 
e-mail: F . VoglOgmx . at 

(Received ) 



Dedicated to Robert F. Tichy on the occasion of his 50 birthday. 



Abstract 

The asymptotic behaviour of the solutions of Poincare's functional equation f(Xz) — 
p(f{ z )) (A > 1) for p a real polynomial of degree > 2 is studied in angular regions W 
of the complex plain. It is known [9j 110] that f(z) ~ exp(z p F(log A z)), if f(z) — > 00 for 
z — > 00 and z £ W, where F denotes a periodic function of period 1 and p = log A deg(p). 
In the present paper we refine this result and derive a full asymptotic expansion. The 
constancy of the periodic function F is characterised in terms of geometric properties of 
the Julia set of p. For real Julia sets we give inequalities for multipliers of Pommerenke- 
Levin-Yoccoz type. The distribution of zeros of / is related to the harmonic measure on 
the Julia set of p. 
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1. Introduction 

1-1. Historical remarks 

In his seminal papers [361 137] H. Poincare has studied the equation 

f(Xz) = R(f(z)), zeC, (1.1) 

where R(z) is a rational function and A 6 C. He proved that, if R{0) = 0, R'(0) = A, and 
|A| > 1, then there exists a meromorphic or entire solution of (|l.lj) . After Poincare, (II. ip 
is called the Poincare equation and solutions of are called the Poincare functions . 
The next important step was made by G. Valiron |45l I46j . who investigated the case, 
where R{z) = p(z) is a polynomial, i.e. 

f(\z)=p(f(z)), zeC, (1.2) 

and obtained conditions for the existence of an entire solution f(z). Furthermore, he 
derived the following asymptotic formula for M{r) — max| z |< r |/(z)|: 

log M(r) - r p F ( ^j- ) , r -> oo. (1.3) 



log|A| 

Here F(z) is a 1-periodic function bounded between two positive constants, p = ^° g |^ 
and d = degp(z). 

Different aspects of the Poincare functions have been studied in the papers [9j 1101 1121 
1141 1211 141] . In particular in [9], in addition to (|1.3[) . asymptotics of entire solutions f{z) 
on various rays re 1 ® of the complex plane have been found. 

It turns out that this asymptotic behaviour heavily depends on the arithmetic nature 
of A. For instance, if arg A = 2tt(3, and j3 is irrational, then f(z) is unbounded along any 
ray argz = i9 (cf. [9]). 

1-2. Assumptions 

In the present paper we concentrate on the simplest, but maybe most important case 
for applications, namely, when A is real and p(z) is a real polynomial (i. e. all coefficients 
of p(z) are real). 

It is known from [46] and [9] that, if f(z) is an entire solution of (|1.2[) . then the only 
admissible values for /o = /(O) are the fixed points of p{z) (i. e. p(fo) = /o). Moreover, 
entire solutions exist, if and only if there exists an % £ N such that 

A"" = P \h). 

It was proved in [9j Propositions 2.1-2.3] that the general case may be reduced to the 
simplest case 

/ (0) = p(0) = and p'(0) = A > 1 

by a change of variables. In the same vein, we can assume without loss of generality that 
/'(0) = 1 and the polynomial p is monic (i. e. the leading coefficient is 1) 

p(z) = z d +p d _ 1 z d - 1 + ---+p 1 z. (1.4) 

1-3. Poincare and Schroder equations 

The functional equation (|1.2[) with the additional (natural) conditions f(0) = and 
/'(0) = 1 is closely related to Schroder's functional equation (cf. [42] ) 

9 {p{z))=\g{z), g(0) = and g'(0) = 1 (1.5) 
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which was used by G. Koenigs [22LI23] to study the behaviour of p under iteration around 
the repelling fixed point z = 0. By definition, g is the local inverse of / around z = 0. 
Both functions together provide a linearisation of p around its repelling fixed point z = 

g(p(f(z))) = Xz and ff(p<*> (/(*))) = X n z, 

where p^ n \z) denotes the n-th iterate of p given by p^ (z) — zandp' n+1 '(z) = p(p^ n \z)). 

We note here that (jTTTJ) and (|1.2p are also called Schroder equation by some authors. 
For instance, the value distribution of solutions of the Poincare (alias Schroder) equation 
has been investigated in |21j . 

1-4. Branching processes and diffusion on fractals 

Iterative functional equations occur in the context of branching processes (cf. |20j). 
Here a probability generating function 

oo 

q{z) = J2 PnZn 

encodes the offspring distribution, where with p n > is the probability that an individual 
has n offspring in the next generation (note that q(l) = 1). The growth rate A = q'(l) 
decides whether the population is increasing (A > 1) or dying out A < 1. In the first case 
the branching process is called super- critical. The probability generating function q^ n \z) 
(n-th iterate of q) encodes the distribution of the size X n of the n-th generation under 
the offspring distribution q. In the case of a super-critical branching process it is known 
that the random variables X~ n X n tend to a limiting random variable Xrx,. The moment 
generating function of this random variable 

f(z) = Ee~ zX ~> 

satisfies the functional equation (cf. [20] ) 

f(Xz) = q(f(z)), 

which is (|1.2[) . if q is a polynomial. Furthermore, this equation can be transformed into 
(|1.2[) . if q is conjugate to a polynomial by a Mobius transformation, especially q{z) = 
^TTz\ ' wnere P is a polynomial. 

Branching processes have been used in [1] [2J 130] to model time for the Brownian 
motion on certain types of self-similar structures such as the Sierpifiski gasket. In this 
context the zeros of the solution of (|1.2p are the eigenvalues of the infinitesimal generator 
of the diffusion ("Laplacian"), if the generating function of the offspring distribution is 
conjugate to a polynomial (cf. |10| I18| I31[ 143 [144]). In this case the zeros of / have to be 
real, since they are eigenvalues of a self-adjoint operator. This motivates the investigation 
of real Julia sets in Section 2J 

1-5. Contents 

The paper is organised as follows. 

In Section |2~T1 we study the asymptotic behaviour of f(z) in those sectors W of the 
complex plane, where 

f(z) — > oo for z — > oo, z eW. (1.6) 
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It was proved in [9j 110) that (|1.6[) implies 

f(z) ~ exp ( z p F ] for z -> oo, z e W, 



log A 

where -F(z) is a periodic function of period 1. In Section l2Tl we will refine this result to 
a full asymptotic expansion of f(z), which takes the form 

oo 

f(z) = exp {z"F (log A z)) + J2 cn exp {-nz"F (log A z)) , (1.7) 

n=0 

where F is a periodic function of period 1 holomorphic in some strip depending on W 
and p = log A d. The proof is based on an application of the Bottcher function at oo of 
p(z). 

We note here that E. Romanenko and A. Sharkovsky [41) have studied equation 
(11.2)) on R (rather than C) and obtained a full asymptotic expansion of this type by 
Sharkovsky's method of "first integrals" or "invariant curves" . 

Further analysis of the periodic function F occurring in 1)1. 7)) is presented in Section l2~3) 
where the Fourier coefficients of F are related to the Bottcher function at oo of p(z) and 
the harmonic measure on the Julia set of p. 

In Section |2~^H the asymptotic behaviour of f(z) is studied in sectors that are related 
to basins of attraction of finite attracting fixed points. 

In Section [3] we relate geometric properties of the Julia set to the location of the zeros 
of/. 

Section Q] is devoted to the special case of real Julia sets J{p). Here we prove, in 
particular, the following inequalities of Pommerenke-Levin-Yoccoz type for multipliers of 
fixed points £: 

= f \p'(0\ > d for mm J(p) < £ < max J(p) 

I > d 2 for £ = minjT'(p) or £ = maxj(p). 

Furthermore, equality can hold only, if p is linearly conjugate to a Chebyshev polynomial 
of the first kind. 

In Section[5]we continue the study of Dirichlct generating functions of zeros of Poincare 
functions that we started in |10j in the context of spectral zeta functions on certain 
fractals. We relate the poles and residues of the zeta function of / to the Mellin transform 
of the harmonic measure fi on the Julia set of p. Furthermore, we show a connection 
between the zero counting function of / and the harmonic measure /i of circles around 
the origin. 



2. Relation of complex asymptotics and the Fatou set 

Throughout the rest of the paper we will use the following notations and assumptions. 
Let p be a real polynomial of degree d as in (|1.4|) . We always assume that p(0) = 
and p'(0) = a\ = A with |A| > 1. We refer to p2 133) as general references for complex 
dynamics. 

We denote the Riemann sphere by C M and consider p as a map on Coo • We recall that 
the Fatou set !F(p) is the set of all z £ Coo which have an open neighbourhood U such 
that the sequence (p^)neN is equicontinuous on U in the chordal metric on Coo- By 
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definition T(p) is open. We will especially need the component of oo of T(p) given by 

Foo{p) = \z G C | lim = 00} , (2.1) 

as well as the basins of attraction of a finite attracting fixed point wq (p(wq) = Kio, 

b'K)l<i) 

T wa (p) = \z e C I lim p<">(*) = too} . (2.2) 

The complement of the Fatou set is the Julia set J(p) — \ T(p). 
The filled Julia set is given by 

K[p) = [z e C I (p {n) (z)) neN is bounded} = C \ ^(p). (2.3) 

Furthermore, it is known that (cf. |15j ) 

0£(p) = dTUp) = J{ P ). (2.4) 

In the case of polynomials this can be used as an equivalent definition of the Julia set. 
We will also use the notations 

W Qi/3 = {zeC\{0}|a<argz</3} (2.5) 

and 

B(z, r) — {w e C I \z - w\ < r} . 

2-1. Asymptotics in the infinite Fatou component 

In O HOj the asymptotics of the solution of the Poincare equation (|1.2p was given. We 
want to present a different approach here, which gives a full asymptotic expansion. 

Theorem 2-1. Let f be the entire solution of the Poincare equation (|1.2p for a real 
polynomial p with A = p'(0) > 1. Assume further that the Fatou component of 00, JFoo(p) 
contains an angular region W a ,p- 

A Then the following asymptotic expansion for f is valid for all z G W a ,p large enough 

00 

f{z) = exp (z p F (log A z)) + J2 c n exp (—nz p F (log A z)) , (2.6) 

71=0 

where F is a periodic function of period 1 holomorphic in the strip 

zeC|A<3z< ' 



log A log A 

and p — log A d. Furthermore, 

Vz G W a ,p : 5Rz p F(log A z) > (2.7) 

holds. 

B Let g denote the Bottcher function associated with p, i. e. 

(g(z)) d = g{p{z)) (2.8) 
in some neighbourhood of 00. Its inverse function is given by the Laurent series around 



g ^H= W + Y,^- (2-9) 
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Then we have 

f{z)=g^ (exp(z'F(log A z))) 
and c n can be determined from the coefficients of p. 

Proof. We recall that p has a super- attracting fixed point of order d = degp at infinity. 
We consider the Bottcher function g associated with this fixed point (cf. O G>J EH 1261 ). 
which satisfies the functional equation (12. 8|) in some neighbourhood of infinity. The 
Bottcher function has a Laurent expansion around infinity given by 

(2 - 10) 

n=0 

which converges for \z\ > R for some R > 0. The coefficients (6 n )neN can be determined 
uniquely from the coefficients of the polynomial p. 

Using the Bottcher function we can rewrite the Poincare equation assuming that 

\m\>R 

(g(f(z))) d = g(p(f(z))) = g(f(\z)). (2.11) 
From this we derive that h(z) = g(f(z)) satisfies the much simpler functional equation 

(h(z)) d = h(\z), 

which only holds for those values z for which \f(z)\ > R. This equation has solutions 

h(z)=exp(z>>F(log x z)) (2.12) 

with p = log A d and F a periodic function of period 1 holomorphic in some strip parallel 
to the real axis. Since \h(z)\ > 1 for all z with \f(z)\ > R by the properties of the function 
g, we have ()2.7|) . 

By (I2.10p g is invcrtible in some neighbourhood of oo and we can write (|2.9[) where the 
coefficients c n depend only on the coefficients of the polynomial p. This function satisfies 
the functional equation 

g(- 1 \w d )=p(g(- 1 Xw)) (2.13) 

for w in some neighbourhood of oo. Inserting (|2.12p into (|2.9|) yields (|2.6p giving an exact 
and asymptotic expression for f{z). □ 

Remark 2-1. E. Romanenko and A. Sharkovsky have studied equation (|1.2|l on K 
(rather than on C) in [41] . Applying Sharkovsky's method of "first integrals" ("invariant 
graphs") they obtained a full asymptotic formula of type (12. 6|) for all solutions f(x), such 
that f(x) — > oo for x — > oo. 

2-2. Bottcher functions, Green functions, and constancy of the periodic function F 
We will make frequent use of the integral representation of the Bottcher function 

g(z) = exp / log(z - x) dp(x) , (2.14) 
\ J J(p) J 

where fj, denotes the harmonic measure on the Julia set J{p) (cf. [H [7j I39j ). This shows 
that g is holomorphic on any simply connected subset of J^ooip)- The measure fj, can be 
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given as the weak limit of the measures 

^» = jr E ( 2 - 15 ) 

p(»)(x)=€ 

where £ can be chosen arbitrarily (not exceptional) and S x denotes the unit point mass 
at x (cf. [71139]'!. 

The function g(z) can be continued to any simply connected subset U of Coo \ fc(p) 
(this follows for instance from the integral representation (|2.14p ). Furthermore, it follows 
from [3j Lemma 9.5.5] and (l2~8]) that 

g(U) C {zeCoo | \z\ > 1}. 

The function log \g(z) \ is the Green function for the logarithmic potential on Foaip) (cf- P2 
Section 9]). Combining classical potential theory with polynomial iteration theory we get 

lim \g{z)\ = 1 ^ z Q E J{p), (2.16) 

Z 'Zq 

where the implication <= is [3j Lemma 9.5.5]. The opposite implication is a general 
property of the Green function (cf. [161 Chapter III], and [391 Section 6.5]) combined 
with the fact that dJ r OQ (p) = Jip) for polynomial p. 

Theorem 2-2. The periodic function F occurring in the asymptotic expression (]2.6p 
for f is constant, if and only if the polynomial p is either linearly conjugate to z d or to 
the Chebyshev polynomial of the first kind Td{z). 

Proof. The periodic function F is constant, if and only if the function h(z) = g(f(z)) 
introduced above satisfies 

h(z) = exp (Cz») (2.17) 

for some constant C ^ 0. This implies that for any wq E J{p) \ {0} the function g has an 
analytic continuation to some open neighbourhood of wq. Thus (|2. 16[) can be replaced 
by ' 

\g(w )\ = l«w e Jip) 

in our case. By (|2. 1T|) this is equivalent to wq = /(zo) for Cz^ S iR. Since Cz p E iM. 
describes an analytic curve (with a possible cusp at z = 0), the Julia set of p is the image 
of this curve under the entire function /, thus itself an analytic arc. 

By |191 Theorem 1] Jip) can only be an analytic arc, if the Julia set of p is either a line 
segment or a circle. The Julia set is a line segment, if and only if p is linearly conjugate 
to the Chebyshev polynomial Td (cf. [3] Theorem 1.4.1]); the Julia set is a circle, if and 
only if p is linearly conjugate to z d (cf. [3j Theorem 1.3.1]). □ 

Remark 2-2. Suppose that the periodic function F is constant. If p is linearly conjugate 
to a monomial, then the Bottcher function g and therefore its inverse are linear functions. 
In this case p = 1. (We recall that we generally assume that /'(0) = 1.) If p is linearly 
conjugate to a Chebyshev polynomial, g^" 1 ^ is linearly conjugate to the Joukowski func- 
tion z H — . In this case p = 1, if is an inner point of the line segment Jip), and p = ^, 
if is an end point of the line segment Jip) (cf. Sections 14-11 and l 4-2[) . Furthermore, the 
asymptotic series (12 .6[) is finite, if the periodic function F is constant. 
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2-3. Further analysis of the periodic function 

In this section we relate the periodic function F occurring in ()2.6|1 to the local behaviour 
of the Bottcher function at the fixed point /(0) = 0. 

This will allow to express the Fourier coefficients of F in terms of residues of the Mellin 
transform (cf. |11U34| ) of the harmonic measure fj, given by (|2.15| . This Mellin transform 
was introduced and studied in [4]. A similar relation was also used in [18] to derive an 
asymptotic expression for / in a special case. 

We will use the relation 

G(w) = log g(w) = / log(> - i) dfi(x) (2.18) 
JJ{P) 

between the (complex) "Green function" G and the Bottcher function g. Assume that 
the Fatou component Too{p) contains an angular region centred at the fixed point 0. 
Furthermore, assume that l\m w ^$ g(w) = 1. Then (|2.12jl holds in this angular region. 
This fact can be used to analyse the local behaviour of \ogg{w) around w — 0: 

log g(w) = ^-^(w)Y F (log x f^(w)j =w"F(log x w)+0(wP +1 ). (2.19) 

Thus the behaviour of the Green function G at the point exhibits the same periodic 
function F as the asymptotic expansion of log / around oo. 




Fig. 1. Paths of integration. 

We now relate the Green function G(w) to the Mellin transform of /x 

M^s) = f {-x) s d»(x), (2.20) 
JJ(p) 

where the branch cut for the function (— x) s is chosen to connect with oo without any 
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further intersection with Following the computations in [4j Section 5] we obtain 

M ^ (s) = hi i [ ' z)s dG{z) = id £ ( ~ z)s dG[z) - 

For 3ts < we have for the circle of radius R 



r (- Z )'dG(z) 

2tt« J\z\=r 



l\z\=R 

which allows to let R — > oo in this case. This gives 



^00 = ^1/ {-z) s dG(z)-j (-zydG(z) 



2m 



2ni 



x s G'(x)dx = s^^ I x s - 1 G(x)dx, 



which relates the Mellin transform of the measure p to the Mellin transform of the 
function G(z) 



MG(s) = I x^Gix) dx = — MJs) for - p < 3?s < 0. (2.21) 

Jo SSUITTS 

The function M M (s) (and therefore M.G(s) by (|2.21[) ) has an analytic continuation by 
the following observation 

M ^ s ) = lT.i (-P^H^yd^x), (2.22) 

where pi ^ (fc = 1, . . . , d) denote the d branches of the inverse function of p; we choose 
the numbering so that p[ (Q) = 0. The summands for k = 2, . . . , d are clearly entire 
functions in s, since the integrand is bounded away from and oo. For the summand 
with k = 1 we observe that 

p[~ 1] {x) = \x + 0{x 2 ) for x -> 0. (2.23) 
A 

Inserting this into (|2.22[) gives 

M^s) = -\- s [ {-x) s dp(x) + -X~ s [ {-x) s O{x) dfi(x) 
a Jj(p) a J J{p ) 



IT,! (-pi'^^Yduix), 



d 

where the second term on the right-hand-side originates from inserting the holomorphic 
function 0(x 2 ) from (I2.23|) into the integrand, which gives a function holomorphic in a 
larger domain. Thus we obtain 

M^s) = — L-^s) (2.24) 

for some function H(s) holomorphic for Sis > — p— 1 (p = log A d). The numerator d\ s — 1 
has zeros at s = —p + (fceZ), which give possible poles for the function M^s). 

Remark 2-3. Using the full Taylor expansion of Pi (x) instead of the O-term in (I2.23|) 
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would yield the existence of a meromorphic continuation of M^{s) to the whole complex 
plane. 

Taking (|2.21|) and (|2.24[) together gives the analytic continuation of M.G(s) to — p— 1 < 
5is < 0. Then the Mellin inversion formula (cf. [11] ) gives (for — p < c < 0) 

c+ioo c+zoo 

TT 1 

' ds 



G(x) = _L /" MG(s)x~ s ds = J^ { 1 - H(s)a 

w 2ni J w 2wi J s sin its d\ s - 1 v ; 

c— ioc c—ioo 

— p— ^-+200 

= — / — H(s) X - s ds + V Res 7WG(s)a;- s . (2.25) 

27ri 7 SSin7TsdA s -l ^ s= _ (3+ 2fcrj 



— P— 77—lOG 



The integral in the second line is 0(x p+ ^), the sum of residues can be evaluated further 
to give the Fourier expansion of the periodic function F 

V Res M G(s)x~ s ^ x"Y f k e 2kmlo ^ x = x p F(\og x x). (2.26) 

The Fourier coefficients fk are given by 

fk = Res . MG(s) = * — T Res . M„(«) 



(-logd-2A7«)sin7r(-p-^) V lo S A 



log A, 

2-4. Asymptotics in a finite Fatou component - analysis of asymptotic values 

It is clear from the functional equation (|1.2p for / that any asymptotic value of / has 
to be an attracting fixed point of the polynomial p (including oo). Thus the analysis in 
Section 1211 can be interpreted as the behaviour of / when approaching the asymptotic 
value oo. In the present section we extend this analysis to all asymptotic values. 

First we study the case of a finite attracting, but not super-attracting fixed point. Let 
Wo be an attracting fixed point of p and denote r\ = p'(wo) (\n\ < 1). Then there 
exists a solution 'J of the Schroder equation 

r?St(z) = y{p{z)), *(wo) = 0, and = 1, (2.28) 

which is holomorphic in !F Wo (p) (for instance, the sequence {rj~ n (jp- n ' (z) — wo))neN con- 
verges to ^ on any compact subset of T Wo {p)). Assume now that T Wo (p) contains an 
angular region W a ,f) H 5(0, r) for some r > 0. Then by conformity of / some angular 
region at the origin is mapped into W a _p n 5(0, r). We consider the function 

j(z) = *(/(*)), 

which satisfies the functional equation 

j(Xz) = *(/(Az)) = *(?(/(*))) = »7*(/(a0) = m{z). (2.29) 

This equation has the solution 

j(z)=z lo ^H(log x z) (2.30) 
with some periodic function of period 1, holomorphic in some strip. This periodic function 
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can never be constant, since otherwise j(z) would have an analytic continuation to the 
slit complex plane. From this it would follow that / is bounded in the slit complex plane, 
a contradiction. 

The function '5 has a holomorphic inverse around 

oo 
n=2 

which allows us to write 

oo 

f(z) = (z l °^ "ff(log A z)) = wo + z l °^ "ff(log A z) + V^" log ^(ff Qog A z)) n , 

n=2 

(2.31) 

which is valid in the angular region W a< p for z large enough. This gives an exact and 
asymptotic expression for / in an angular region. 

In the case of a super-attracting fixed point u>o we have p'(wo) — 0. Assume that 
the first k— 1 derivatives of p vanish in wq, but the fc-th derivative is non-zero. Then 
p(z) = (z — wo) k P(z) with P(wq) = A ^ 0. We use the solution g of the corresponding 
Bdttchcr equation 

g(p(z)) = A(g(z)) k g(w )=0, g'(w Q ) = l (2.32) 

to linearise (11.21) 

g(f(\z))=g(p{f(z)))=A(g{f(z))) k . 
Thus the function h(z) = g(f(z)) satisfies 

h(\z) = A{h(z)) k . 

This equation has solutions 

h(z) = A-t^ cxp (z log * k L (log A z)) 

for a periodic function L of period 1 and a suitable choice of the (k — l)-th root. Fur- 
thermore, by the fact that lim^oo h(z) = we have 

R (z l °^ k L (log A zj) < for f(z) e F m (p). 

using the local inverse of g around we get 

f(z) = g(-V (a~& exp (z'°^ k L (log A *))) 

= w Q + A-t^ exp (z log * k L (log A z)) (1 + o(l)). (2.33) 

Summing up, we have proved 

Theorem 2-3. Letwo be an attracting fixed point ofp such that the Fatou component 
T Wo (j>) contains an angular region W ay p fl 5(0, r) for some r > 0. Then the asymptotic 
behaviour of f for z — > oo and z G W a ,p is given by (|2 . 3 1 1) . if r] = p'(wo) ^ 0, and by 
(|2.33p . if p(z) — wq has a zero of order k in Wq. 

Remark 2-4. The periodic function H in (|2.31[) cannot be constant, because otherwise 
f(z) would be bounded. The periodic function L in (|2.33j) can only be constant, if p is 
linearly conjugate to z k , by the same arguments as in the proof of Theorem l2-2l (the case 
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of Chebyshev polynomials does not occur, because they only have repelling finite fixed 
points). 

As a consequence of Ahlfors' theorem on asymptotic values (cf. [17p and Valiron's 
theorem on the growth of / (cf. [451 [46]) we get an upper bound for the number of 
attracting fixed points of a polynomial. 

Theorem 2-4. Let p be a real polynomial of degree d > 1 and let 

7 = max{\p'(z)\ | p(z) = z) . 

Then the number of (finite) attracting fixed points of p is bounded by 21og 7 d, i.e. 

# {z e C | p(z) = z A \p'(z)\ < 1} < 2 log 7 d. (2.34) 

3. Zeros of the Poincare function and Julia sets 
In this section we relate the distribution of zeros of the Poincare function in angular 
regions to geometric properties of the Julia set 3{p) of the polynomial p. 

Theorem 3T. Let p be a real polynomial with p(0) = and p'(0) = A > 1. Then the 
following are equivalent 

(i) Vr > : W a ,p n J(p) n B(Q, r) £ 

(ii) W a ,0 contains a zero of f. 

(iii) W a fi contains infinitely many zeros of f . 



Proof. We first remark that (ii) and (iii) are trivially equivalent, since /(zo) = implies 
that f(X n z ) = 0. 

For the proof of ' (i) (ii) ' we take < e < so small that 

Vr > : W a+£ ,/3- e n J(p) n B(0, r) ± 0. 

Then we take r > so small that 

iy a+E , N nB(o,r)c/(M, (3.1) 

which is possible by conformity of / and /'(0) = 1. Since the preimages of are dense 
in J(p), there exists r] S W a +s,f}-e H 5(0, r) and n G N such that p {n \rf} = 0. By (|3~T|) 
there exists ^ € W a< p such that /(^) = rj, from which we obtain 

/(A n 0=P (n) (/(0)=P (n) fa)=°- 



For the proof of ' (iii) =>■ (i) ' we take z £ W a ^ with f(z Q ) = 0. Then 



For any r > and n large enough f(X~ n zo) G W Q ,/3 n B(0, r), which gives (i) □ 
Similar arguments show 

Theorem 3-2. Let p be a real polynomial with p(0) — and p'(0) = A > 1. Then 

3{p>) C M _ U {0} <f> all zeros of f are non-positive real (3-2) 

and 

3{p) Clo all zeros of f are real. (3-3) 



Poincare functions and Julia sets 13 
4. Real Julia set 

Lemma 4-1. Let p be a real polynomial of degree d > 1. Then the Julia-set J(p) is 
real, if and only if there exists an interval [a, b] such that 

pt-V([a,b])C[a,b]. (4.1) 

Proof. Assume first that J{p) C R and take the interval [a, b] = [min J(p), maxj"(p)] . 
Let e > 0. Since J(p) is perfect, there exist £, 77 G J{p) with a < £ < a+e < fe— e < 77 < 6. 
All preimages of £ and 77 are in J^p) by the invariance of J(p). Furthermore, all these 
preimages are distinct. Therefore, every value x G [£, 77] has exactly d distinct preimages 
in [a, b] by continuity of p. Since e was arbitrary and the two points a, b also have all 
their preimages in J(p) C [a, b], we have proved (14. 1|) . 

Assume on the other hand that [a, b] satisfies (|4.1|) . Since the map p has only finitely 
many critical values, there exists x G [a, b] such that the backward iterates of x are dense 
in the Julia set. By (|4.1|1 all these backward iterates are real; therefore J{p) is real. □ 

Remark 4-1. By the above proof we can always assume [a, 6] = [min (p), max J^(p)]. 
Furthermore, we have 

p ({mm J(p), max J(p)}) C {min,7(p),niax,7(p)}, 

which implies that at least one of the two end points of this interval is either a fixed 
point, or they form a cycle of length 2. 

Theorem 4-1. Letp be a polynomial of degree d > 1 with real Julia setj(p). Then for 
any fixed point £ of p with mm J'(p) < £ < max J(p) we have ^ d- Furthermore, 

\p'(mmj(p))\ > d 2 and \p' (max J' (p))\ > d 2 . Equality in one of these inequalities implies 
that p is linearly conjugate to the Chebyshev polynomial Td of degree d. 

Remark 4-2. This theorem can be compared to [8, Theorem 2] and [291 138) . where 
estimates for the derivative of p for connected Julia sets are derived. Furthermore, in 
[13] estimates for — log \ (p^)'(z)\ for periodic points of period n are given. 

Before we give a proof of the theorem, we present a lemma, which is of some interest on 
its own. A similar result is given in |27l Chapter V, Section 2, Lemma 3]. 

Lemma 4-2. Let f be holomorphic in the angular region W a ,p If there exists a positive 
constant M such that 

Vz&W a , : \f(z)\>M, 

then 

V£> 3A,B > Vz G W a+e ^ £ : \f(z)\ < Bexp(^H 

with K — tt—- 

p—ot 

Proof. Without loss of generality we can assume that M — 1 , a = — ^ , and /3 = ^ . In 
this case k = 1. The function 

v(z) = \og\f(z)\ 

is a positive harmonic function in the right half-plane. Thus it can be represented by the 
Nevanlinna formula (cf. [281 p. 100]) 

, x f°° dv(t) 
v(x + iy) = - y -^ + <rx, (4.2) 
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where v denotes a measure satisfying 

du{t) 



1 + t 2 



< oo 



and a > 0. 

In the region given by | argz| < 5 — £ and \z\ > 1 we have 

\z — it\ > max(|i| sine, |z| sine) > max(l, \t\) sine. 
From this it follows that 



which gives 



it\ 2 > -(1 + i 2 )sin 2 e, 



di/(t) < _2_ [°° Mti < B 



it\ 2 ~ sm 2 ej^ 1 + t 2 



for \z\ > 1 and some B £ > 0. Setting A — ^B e + a and observing that x < \z\ completes 
the proof. □ 



Proof of Theorem \4-l\ Without loss of generality we may assume that the fixed point 
£ = 0. Then we consider the solution / of the Poincare equation 

f(\z)=p(f(z)) 

with A = p'(O). We assume first that A > 0. 

First we consider the case minj'(p) < £ < vaaxj(p). In this case the function f{z)/z 
tends to infinity uniformly for z — > oo in the region e < arg z < n — e for any e > by 
Theorem 12- II Furthermore, we know that 

\f(z)\ > Cexp(A\z\ lo ^ d ) 

in this region for some positive constants A and C. Since f(z)jz does not vanish at z = 0, 
this function satisfies the hypothesis of Lemma 14-21 from which we derive that 

n 



log A d < 



it - 2e 

holds for any e > 0, which implies A = p'(0) > d. 

The proof in the case £ = maxj7Xp) runs along the same lines. The function f(z)/z 
tends to infinity uniformly in any region | argz| < 7r — e in this case, which by Lemma 14-21 
implies 

i 7 "ft 
log, d < 

for all e > 0, and consequently A = p'(Q) > d 2 . 

For negative A = p'(0) we apply the same arguments to p^ 2 \ 

For the proof of the second assertion of the theorem, we first assume that the fixed 
point £ = satisfies a = mvaj(p) < < maxj7(p) = b and that p'(0) = d. We know 
that for a suitable linear conjugate q of the Chebyshev polynomial Td we have q'(0) = d 
and J(q) = [a, b] with G (a, b). 

Let us assume now that p'(0) = d and J(p) is a Cantor subset of the real line, or 
after a rotation that J{p) is a Cantor subset of the imaginary axis (this makes notation 
slightly simpler). 
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By arguments, similar to those in the beginning of Section [2-31 we can write 

H(z)=mogg(f(z))= [ log\f(z)-x\dti(x). (4.3) 
Jj(p) 

Since 5Rlogg(.) is the Green function of J(p) with pole at oo (cf. [3l Lemma 9.5.5] or 
[39]). we know that H(z) > for all z G C and H(z) = 0, if and only if f(z) G J(p) 
(since IC(p) = J(p) in the present case). By Theorem 12 -II we have 

H(z) = K (zF(log d z)) = xK{F{\og d z)) - y%{F{\og d z)) for z = x + iy, (4.4) 

and by Theorem 12-21 the function F is not constant in the present case. The periodic 
function $sF(t + if) has zero mean, since the mean of F is real. Thus QF(t + itp) attains 
positive and negative values for any ip. We now take z = iy G zR + to obtain 

H(iy) = -yS(F(log d y + i^^)). 

Since SsF attains positive values by the above argument, we get a contradiction to H(z) > 
for all z. 

A similar argument shows that for = maxj'(p) and p'(0) — d 2 the assumption that 
the Julia set is not an interval leads to the same contradiction. □ 

Remark 4-3. Lemma 6.4 in [9] proves Theorem 14 II for the special case of quadratic 
polynomials. The proof given in [9] is purely geometrical. 

Remark 4- 4. We have a purely real analytic proof for \p' (max J{p))\ > d 2 , which is 
motivated by the proof of the extremality of the Chebyshev polynomials of the first kind 
given in 40 1. However, we could not find a similar proof for the other assertions of the 
theorem. 

4-1. The Julia set is a subset of the negative reals 

As a consequence of Lemma 14-21 we get that any solution of the Poincare equation for 
a polynomial with Julia set contained in the negative real axis has order < h. The only 
solutions of a Poincare equation with order ^ in this situation are the functions 

f(z) = — ^cosh \j2az — lj 

for 

p(z) = (T d {az + 1) - l)/ffl, 

where a G IR + and T<i denotes the Chebyshev polynomial of the first kind of degree 
d. This is also the only case where the periodic function F in (|2.6p is constant in this 
situation. 

Corollary 4-1. Assume that p is a real polynomial such that S{p) is real and all 
coefficients pi (i > 2) of p are non-negative. Then J{p) C R - U {0} and therefore 

/( 2 )~exp(^(g£)) (4.5) 

for z — > oo and \ argz| < tt. Here F is a periodic function of period 1 holomorphic in the 
strip given by \ < 5sw\ < j^j- Furthermore, for every e > Sfte' pa '' sz J 7 ' ( }°g^ ) is bounded 
between two positive constants for \ argz| < tt — e. 
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Proof. From [9l Lemmas 6.4 and 6.5] it follows that f(z) has only non-positive real 
zeros. Then by Thcorem l3-2l .7 (v) C R~ U {0}. Finally, the assertion follows by applying 
[9j Theorem 7.5]. □ 

Example 1. In order to illustrate the above results, we shall turn to the equation 

/(5z) = 4/(z) 2 -3/(z), 

which arises in the description of Brownian motion on the Sicrpihski gasket [101 1241 1251 
143] , Here p(z) — 4z 2 — 3z, and the fixed point of interest is /(0) = 1. This fits into the 
assumptions of Section [L2l only after substituting g(z) — 4(/(z) — 1), where g satisfies 

g(5z) = g{zf + 5g{z). 

Now Corollary 14-11 may be applied to this equation (the preimages of are real by [9] 
Lemma 6.7]) to give (|4.5j) . 

Note also that p'(Q) — 5 > 4 = 2 2 in accordance with Theorem 14- II 

4-2. The Julia set has positive and negative elements 

Again as a consequence of Theorem 14-11 the solution of the Poincare equation for a 
polynomial with real Julia set with positive and negative elements has order < 1. The 
only solution of a Poincare equation of order 1 in this situation are the functions 

for 

p(z) = - (T d (a(z + &)) - 60 , 
a 

where a £ R + and = cos jrj f° r 1 < fc < ■ This is again the only case where the 
periodic function F in (|2.6| is constant in this situation. 

5. The Zeta function of the Poincare function 

In |10j the zeta function of a fractal Laplace operator was related to the zeta function of 
certain Poincare functions. Asymptotic expansions for the Poincare functions were then 
used to give a meromorphic continuation of these zeta functions as well as information on 
the location of their poles and values of residues. In this section we give a generalisation 
of these results to polynomials whose Fatou set contains an angular region W- a , a around 
the positive real axis. In this case the solution / of (|1.2p has no zeros in an angular region 
W- a . a . Furthermore, from the Hadamard factorisation theorem we get 

/« - ^ (h-^'-fj /( n o (i +!)«*(-!+£+•■■+ c-i>"£) . 

(5-1) 

where k — [log^ d\ . By the discussion in 10, Section 5] the values e±, . . . , are given 
by the first k terms of the Taylor series of log 

\o g f -^=Y,(-ir iej f+°(z k+1 )- 
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The zeta function of / is now defined as 

000= E £" s ' ( 5 - 2 ) 

where is dehned using the principal value of the logarithm, which is sensible, since 
£ is never negative real by our assumption on !Foa{p)- The function C/( s ) is holomorphic 
in the half plane Sfts > p. In [10] we used the equation 



k 



log/Or) - log* - Y2(-iy-^ a;— 1 dx = O(s)— — , (5.3) 



=i 



s sin 7rs 



which holds for p < -fts < k + 1, to derive the existence of a meromorphic continuation 
of to the whole complex plane. There ( |10l Theorem 8]) we obtained 

fk ( 2nik\ ( lixik 

Re H . C/(s) = -— P + t— r smvr p 1 



where is given by (|2.27|) . From this we get 

Res . Cf(s) = - Res . M^s). (5.4) 

This shows that the function 

Cf(s) - M^-s) (5.5) 

is holomorphic in p — 1 < 3?s < p + 1, since the single poles on the line 3ts = p cancel. 
This fact was used in |18j to derive an analytic continuation for £/( s )- 

Theorem 5-1. Let f be the entire solution of (|1.2p and assume that p is neither 
linearly conjugate to a Chebyshev polynomial nor to a monomial and that W- a ,a C fooip) 
for some a > 0. Then the following assertions hold 

(i) the limit lim^oo t~ p log f(t) does not exist. 

(ii) C/( s ) has at least two non-real poles in the set p + l-KiaTL [a — lo g ^ ) ■ 
(hi) the limit lim^^o x~ p G(x) with G given by (|2.18|) does not exist. 

Proof. Equation (|2.6p in Theorem 12-11 (see also [9]) implies that 

z _p log/(z) = F(log A z) + o(l) for z — > oo and z £ W- a ,a 

with a periodic function F of period 1 . Theorem 12-21 implies that F is a non-constant . 
Thus the limit in (i) does not exist. 

Since the periodic function F is non-constant, there exists a io / such that the 
Fourier-coefficients f±k do not vanish. By (|2.6j) we have 

log f(z) ^z P Y,f^ m + 0(z- M ) 
feez 

for any M > 0. By properties of the Mellin transform (cf. |35j ). every term Az p+lT in the 
asymptotic expansion of log f(z) corresponds to a first order pole of the Mellin transform 
of log/(z) with residue A at s = p + it. Since fk ^ 0, from (|5.3[) we have simple poles 
ofC/(s)at S = p±f^ 



Assertion (iii) follows from (i) by (|2.19[) . □ 
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In the following we consider the zero counting function of / 



N f (x)= ]T 1. (5.6) 



\£\<x 
/(«)=0 



Theorem 5-2. Let f be the entire solution of (|1.2[) . Then the following are equivalent 

(i) the limit lhn x ^ toc x~ p Nf(x) does not exist. 

(ii) the limit \hn t ^ot~ p fi(B(0,t)) does not exist. 



Proof. For the proof of the equivalence of (i) and (ii) we observe that by the fact that 
/'(0) = 1, there is an ro > such that / : B(0, ro) — > C is invertible. For the following 
we choose n — |log A (:r/ro)J + k and let the integer k > be fixed for the moment. Then 
we use the functional equation for / to get 

N f (x) = #{? | /(A"£) = p (n) (/(0) = A ICI < x\- n ] =#(p(-»>(0) n /(B(0, xX~ n ))) . 
This last expression can now be written in terms of the discrete measure [i n given in 

N f (x)=d n (i n (f(B(0,x\- n ))). 

By the weak convergence of the measures /i„ (cf. [7\) we get for x — > oo (equivalently 

n — » cxd) 

JV>(s) = d n /x(/(.B(0,a;A-™))) + o(d") = a^(xA- n )-"A*(/(B(0,arA- n ))) + o{x p ). 

By our choice of n we have roA _fe_1 < x\~ n < ro\~ k ~ 1 , which makes the first term 
dominant. From this it is clear that the existence of the limit 

lim x- p Nf(x) = C 

x — >oo 

is equivalent to 

p(f(B(fl,t))) = Ct p for r \- k < t < rnX-^l 
Since k was arbitrary this implies 

At(/(-B(0, *))) = Ct p for < t < r . (5.7) 
It follows from /'(0) = 1 that 

Ve > : 3d > : Vi < 6 : B(0, (1 - e)t) C f(B(0, t)) C 5(0, (1 + e)t). (5.8) 



Thus the existence of the limit in assertion (ii) is equivalent to 

\unt- p ^(f(B(0,t)))=C. 



Thus (i) and (ii) arc equivalent. □ 

Remark 5-1. If J{p) is real and disconnected then the limits in Theorem 15-21 do not 
exist. Furthermore, it is known that the limit 

\imt- p ^f(B(w,t))) = C 

does not exist for /^-almost all w £ J{p) (cf. [321 Theorem 14.10]), if p is not an integer. 

This motivates the following conjecture. 



Conjecture. The limits in Theorem \5-S\ exist, if and only if p is either linearly con- 
jugate to a Chebyshev polynomial or a monomial. 
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